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ABSTRACT 

Introduction: There are now over 30 confirmed loci 
predisposing to rheumatoid arthritis (RA). Studies have 
been largely undertaken in patients with anticyclic 
citrullinated peptide (anti-CCP) positive RA, and some 
genetic associations appear stronger in this subgroup 
than in anti-CCP negative disease, although few studies 
have had adequate power to address the question. The 
authors therefore investigated confirmed RA susceptibility 
loci in a large cohort of anti-CCP negative RA subjects. 
Methods: RA patients and controls, with serological and 
genetic data, were available from UK Caucasian patients 
(n=4068 anti-CCP positive, 2040 anti-CCP negative RA) 
and 13,009 healthy controls. HLA-DRB1 genotypes and 
36 single nucleotide polymorphisms were tested for 
association between controls and anti-CCP positive or 
negative RA. 

Results: The shared epitope (SE) showed a strong 
association with anti-CCP positive and negative RA, 
although the effect size was significantly lower in the 
latter (effect size ratio=3.18, p<1.0E-96). A non-intronic 
marker at TNFAIP3, GIN1/C5orf30, STAT4, ANKRD55/ 
IL6ST BLK and PTPN22 showed association with RA 
susceptibility, irrespective of the serological status, the 
latter three markers remaining significantly associated 
with anti-CCP negative RA, after correction for multiple 
testing. No significant association with anti-CCP negative 
RA was detected for other markers (eg, AFF3, CD28, 
intronic marker at TNFAIP3), though the study power for 
those markers was over 80%. 
Discussion: In the largest sample size studied to date, 
the authors have shown that the strength of association, 
the effect size and the number of known RA susceptibility 
loci associated with disease is different in the two 
disease serotypes, confirming the hypothesis that they 
might be two genetically different subsets. 



INTRODUCTION 

Based on the presence or absence of anticyclic cit- 
rullinated peptide (anti-CCP) antibodies, rheuma- 
toid arthritis (RA) can be classified into anti-CCP 
positive and anti-CCP negative RA. Anti-CCP 
antibodies have been widely shown to be strong 
predictors of disease severity and radiological dam- 
age.^ ^ It is currently a matter of some debate as to 
whether anti-CCP positive and anti-CCP negative 
RA are two distinct entities or represent two differ- 
ent subsets of one and the same disease. Linkage 
and association analysis revealed the shared 



epitope (SE) to be associated only with anti-CCP 
positive RA and not with anti-CCP negative RA.* 
A study in twin pairs has shown that the estimated 
heritability of anti-CCP negative RA is 66% (95% 
CI 21% to 82%), similar to the heritability of anti- 
CCP positive RA, estimated at 68% (95% CI 55% 
to 79%)7 In the same study, the SE was found to 
explain 18% of the genetic component of RA sus- 
ceptibility in anti-CCP positive RA but only 2.4% 
in anti-CCP negative RA. 

Several studies have investigated putative associ- 
ations between different human leucocyte antigen 
(HLA) alleles or single nucleotide polymorphisms 
(SNPs) within the HLA region and predisposition 
to anti-CCP negative RA,^*^"^^ with contradictory 
results. A large meta-analysis, across four European 
populations, found only a weakly significant asso- 
ciation between several HLA alleles and anti-CCP 
negative RA, but observed marked geographical 
differences.^^ Lack of consistency between studies 
might therefore be explained by a different sample 
size and power, by geographical differences in allele 
frequencies and association patterns or by differ- 
ent study designs or definitions of HLA genotypes 
(two vs four-digit typing, different classifications 
for the SE). A large study, performed on Caucasians 
of Northern European descent, investigated several 
HLA-DRBl susceptibility and protective models 
and SE subgroups for association with RA after 
stratification by autoantibody status. Significant 
associations between several HLA-DRBl alleles 
and anti-CCP negative RA were found. Together, 
these findings strengthen the hypothesis that 
genetic factors predisposing to anti-CCP positive 
RA are different to those predisposing to anti-CCP 
negative RA. 

A recent large meta-analysis brought the num- 
ber of confirmed non-HLA RA susceptibility loci to 
31 However, most studies on the identification of 
RA susceptibility loci published to date have been 
performed in largely anti-CCP positive RA cohorts, 
thereby biasing the search of RA susceptibility loci 
towards genetic variants predisposing to anti-CCP 
positive RA. Importantly, most studies have been 
underpowered to identify anti-CCP negative RA 
predictors. 

Very few studies have systematically com- 
pared the genetic basis of anti-citrullinated protein 
autoantibody (ACPA)-positive to ACPA-negative 
RA outside the HLA region. In a genome-wide 
association study (GWAS) recently performed 
in 774 anti-CCP negative RA patients, 1147 anti- 
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CCP positive RA patients and 1079 common controls/^ no SNP 
achieved genome-wide significance in the comparison between 
anti-CCP negative RA and controls, while the PTPNZZ gene was 
associated with anti-CCP positive RA, together with hundreds 
of SNPs located within the HLA locus on chromosome 6. 

Candidate gene association studies and a subsequent meta- 
analysis have confirmed an association of STAT4 polymor- 
phisms with anti-CCP positive and negative RA.'^ IRF5 
polymorphisms have been shown in independent studies to 
be more strongly associated with anti-CCP negative RA than 
with anti-CCP positive RA in Caucasians, while this differ- 
ential association is controversial in Asians.^" No other SNPs 
have been convincingly associated with anti-CCP negative RA. 
Interestingly, the association of PTPNZZ polymorphisms with 
anti-CCP negative RA is controversial, with some investiga- 
tors reporting association in anti-CCP negative patients. For 
example, a recent study investigating the usefulness of data 
derived from electronic health records tested multiple non-HLA 
RA susceptibility markers in 871 anti-CCP positive RA patients, 
378 anti-CCP negative RA patients and 1212 common healthy 
controls.^'' Only PTPNZZ showed an association with anti-CCP 
negative RA with a p value<0.05. Due to the small sample size 
of the anti-CCP negative subgroup, conclusions could be made 
only for SNPs in aggregate rather than for individual SNPs. The 
authors conclude that there is a partial overlap between the 
genetic basis of anti-CCP positive and anti-CCP negative RA. 

We hypothesised that currently known RA susceptibility 
SNPs would show a differential association pattern in anti-CCP 
negative RA compared with anti-CCP positive RA. Therefore, 
we tested the 31 RA confirmed susceptibility loci for associa- 
tion with RA in a dataset comprising between 1935 and 3827 
anti-CCP positive RA patients, between 808 and 1918 anti-CCP 
negative RA patients and between 11468 and 12392 healthy 
controls per genetic marker, the largest number of anti-CCP 
negative patients studied so far in this setting. 

METHODS 

Data collection, SNP selection and genotyping 

Confirmed RA susceptibility SNPs were selected from the large 
meta-analysis by Stahl etal}'^ Most of the SNPs had been already 
genotyped in RA cases and controls as part of several other proj- 
ects conducted at our laboratory. Cohorts and patient charac- 
teristics are presented in the online supplementary table SI; 
genotyping and quality control procedures have been described 
elsewhere.^^ All patients and controls were Caucasians 
originating from the UK satisfying the 1987 American College 
of Rheumatology classification criteria for RA. Only RA cases 
with an available anti-CCP status, as determined with the 
second generation CCP (CCP2) assay, were included in the 
analysis. Several studies/cohorts did not contain any controls. 
Therefore, anti-CCP positive RA cases or anti-CCP negative RA 
cases from different studies were pooled together and compared 
with controls. If a SNP reported in the study of Stahl was not 
available in a cohort, a SNP proxy was selected to fulfil the fol- 
lowing requirements: linkage disequilibrium (r^) >0.90 with the 
original SNP and a maximum of three different SNPs for one 
specific locus across all datasets. If these requirements could not 
be fulfilled for one dataset, it was excluded from the analysis 
for that particular locus. Therefore, different numbers of cases 
and controls were available for analysis for different loci. The 
following HLA-DRBl genotypes were considered as SE alleles: 
'0101' '0102' '0104' '0105' '0107' '0108' '0110' '0111' '040T 
'0404' '0405' '0408' '0409' '0410' '0413' '0416' '0419' '042T 



'0423' '0426' '0428' '0429' '0430' '0433' '0434' '0435' '0438' 

'0440' '0442' '0443' '0445' '1001' '1113' '1126' '1134' '1402' 

'1409' '1413' '1417' '1419' '1420' '1421' '1429' '1430' '143T 
'1432' '1434' '1441' '1446' '1447' '1448'. 

Statistical analysis 

Multinomial logistic regression was applied to compute OR, 
95% CI and p values for association between the minor allele 
at every locus and either anti-CCP-positive (OR^j.p.pjj5;|.;^J or 
anti-CCP-negative RA (OR^. j^gg^jj^g), assuming additivity on 
the log-odds scale (ie, every locus was coded as 0,1 or 2 cor- 
responding to the copy number of the minor allele). The minor 
allele was defined according to the allele frequency in the total 
population, including cases and controls. To test for differences 
between 0\ ^.^^^ and OR^^p.^^^^^^^^, the linear combination 

- P., where p^ is log(OR,,p.p„,,„J and p. is log(OR^^p.„ ) 
was calculated, along with its standard error. This enables a 
p value for the difference in association to be calculated. 

Due to a high pretest probability for association in anti-CCP 
positive RA, p values were not corrected for multiple testing in 
this subset. Based on previously published data, 32 of the 35 
non-HLA SNPs tested have been confirmed to represent inde- 
pendent effects. Therefore, without assuming any prior proba- 
bifity for association in the anti-CCP negative RA subset, which 
is a stringent assumption, the Bonferroni-corrected significance 
threshold would be 1.6E-03. Statistical analysis was performed 
with Stata V.10.1 (Stata Statistical Software: Release 10; Stata 
Corp., College Station, Texas, USA). 

RESULTS 

Thirty-one confirmed RA susceptibility loci,^'* some of which 
contain independent effects, were considered for analysis. 
Together, this represents a total of 36 markers, plus the SE. 
Actual SNPs and proxys are presented for each locus in table 1. 
The combination of datasets available led to a total number of 
4068 anti-CCP-positive, 2040 anti-CCP-negative RA and 13009 
healthy UK controls. The actual number of cases and controls 
varies for every locus from 1935 to 3827 for anti-CCP positive 
RA, from 808 to 1918 for anti-CCP negative RA and from 11468 
to 12392 for healthy controls. Basic cohort characteristics are 
presented in the online supplementary table SI. The results of 
the association analysis are presented in table 2. As expected, 
every locus showed an association with anti-CCP positive RA 
with a p value<0.05. The SE and the corresponding tag SNP 
were highly associated with anti-CCP positive RA with an OR 
of 4.08 and 2.68, respectively. Three other markers at AFF3, 
CDZ8, PTPNZZ and two at the TNFAIP3 locus reached genome- 
wide significance (<5E-08) in anti-CCP positive RA with OR of 
1.17, 1.18, 1.91, 1.29 and 1.45, respectively. By contrast, only 
six non-HLA loci in total reach a p value below 0.05 in anti- 
CCP negative RA: TNFA1P3, GINl/C5orf30, STAT4, ANKRD55/ 
IL6ST, BLK and PTPNZZ. The three last loci remained significant 
at the Bonferroni corrected threshold of 1.6E-03. The SE, and its 
tag SNP, show an association with anti-CCP negative RA with 
an OR of 1.28 and 1.15, respectively. The online supplementary 
table S2 shows association results stratified by rheumatoid fac- 
tor and anti-CCP positivity. The probability of obtaining, by 
chance, at least 7/36 associated loci at a significance level of 0.05 
is 1.8E-03. Therefore, undera prior hypothesis of no association, 
there is still a significant accumulation of RA susceptibility loci 
associated with anti-CCP negative RA. 

Due to the different number of cases in anti-CCP positive and 
negative RA, simply comparing p values would be misleading. 
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Table 1 Single nucleotide polymorphism (SNP) markers and their proxys for the independent 
rheumatoid arthritis susceptibility loci considered for analysis. The shared epitope is not defined by 
SNP markers, but by a list of 4-digit HLA-DRBl alleles described in the methods section 



Chromosome 


Locus name 


Single nucleotide polymorphism markers 


1 


CD2/CD58 


rsl 1586238 




1 


FCGR2A 


rsl 274661 3 


rsl 0494360 


1 


MMELl/TNFRSFU 


rsl 091 0099 


rs3890745 


1 


PTPN22 


rs2476601 


rs6679677 


1 


PTPRC 


rsl 091 9563 


rsl 932435 


2 


AFF3 locus 1 


rsl 160542 


rsl 1676922 


2 


AFF3 locus 2 


rsl 0865035 




2 


CD28 


rsl 980422 




2 


CTLA4 


rs3087243 


rs231804 


2 


REL 


rsl3031237 




2 


SPRED2 


rs934734 


rsl7534670 


2 


STAT4 


rs7574865 




3 


DNASE1L3/PXK 


rsl3315591 


rs9813011 


4 


IL2/IL21 


rs6822844 


rsl3151961 


4 


RBPJ 


rs874040 


rsl 051 7086 


5 


ANKRD55/IL6ST 


rs6859219 




5 


GIN1/C5orf30 


rs26232 


rs35797 


6 


CCR6 


rs3093023 


rs6907666 


6 


HLA-DRBl 0401 tag 


rs6910071 




6 


PRDMl 


rs548234 




6 


TAGAP 


rs394581 


rsl 69858 


6 


TNFAIP3 locus 1 


rs6920220 


rs2327832 


6 


TNFAIP3 locus 2 


rsl 3207033 


rsl 04991 94 


6 


TNFAIP3 locus 3 


rs5029937 


rs5029939 


7 


IRF5 


rsl 0488631 


rsl2531711 


8 


BLK 


rs2736340 




9 


CCL21 locus 1 


rs951005 




9 


CCL21 locus 2 


rs2812378 


rsl 081 41 38 


9 


TRAF1/C5 


rs3761847 




10 


IL2RA locus 1 


rs706778 


rsl 0795791 


10 


IL2RA locus 2 


rs21 04286 




10 


PRKCQ 


rs4750316 


rsl 0796045 


11 


RAG1ARAF6 


rs540386 


rs5030437 


12 


KIF5A/PIP4K2C 


rsl 678542 


rsl 11 72254 


20 


CD40 


rs4810485 


rsl 569723 


22 


IL2RB 


rs3218253 


rs3218258 



rs9653442 



rs3093024 



rs7072793 



rsl 046864 



Therefore, the significance of the observed difference in associa- 
tion between anti-CCP positive and negative RA was addressed 
by computing a comparison OR (or effect size OR) and its 
corresponding p value (column 'comparison' in table 2). This 
allows the classification of the loci into three distinct categories 
(table 3 and figure 1): the HLA-DRBl SE, PTPN22 and one marker 
at TNFAIP3 are associated with anti-CCP positive and negative 
RA, but show a clear differential association with an effect size 
significantly higher in anti-CCP positive RA (category 1). Other 
loci, like C^orfSO or STAT4 are associated with RA irrespective 
of the serological status, with the effect size not differing signifi- 
cantly between subsets (category 2: anti-CCP independent asso- 
ciations). A third category comprises anti-CCP positive specific 
loci with no significant association detected in anti-CCP negative 
RA; however, the p value for the effect size ratio is below 0.05. 
The remaining loci could not be classified into one of these three 
categories, because although they are associated with anti-CCP 
positive RA with a p value<0.05, the effect in anti-CCP-negative 
RA is not significantly different from that in anti-CCP-positive 
RA, nor significantly different from the null. This last situation 
can only be explained by a lack of power. Indeed, for the vast 
majority of markers in category 3, the study has, in the anti-CCP 



negative subgroup, a power between 60% and 90% at the 0.05 
significance level to detect an association of the same effect size 
as observed in the anti-CCP positive subgroup, while power 
drops to between 15% and 50% for most of the unclassifiable 
markers. As an example in category 3, an association of an effect 
size of 1.17 would be detected [oi AFF3 locus i with a power of 
89.3% at the 0.05 significance level for a minor allele frequency 
(MAP) of 45.5% in controls. An association of a larger effect size 
of 1.45 would be detected for the intronic TNFAIP3 marker with 
a power of 81.6% at the 0.05 significance level for a MAP of 
3.5% in controls. The detection power drops for smaller effects 
around 1.13, but is still 64.4% CCLZi locus 2 (bAAV in controls 
34.2%). However, although CDZ/CD58 has been genotyped 
in 1918 anti-CCP negative cases, it is not possible to classify 
it as being associated with anti-CCP positive, negative or both 
because the power to detect an effect of 1.11 in anti-CCP nega- 
tive with a MAP of 24% in controls is only 46.1%. 

Of note, is the fact, that the effect size of all but two loci, 
ANKRD55 and BLK (figure 1), is larger in anti-CCP positive RA 
than in anti-CCP negative RA. However, due to relatively wide 
CIs, the effect size ratio is not significant, so these two loci are 
both classified in category 2. 
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Interestingly, the three known independent effects at TNFAIP3 
are classified in at least two different categories (table 3). At this 
locus, the most profound discordance between anti-CCP posi- 
tive and negative RA is seen for the intronic marker (figure 1), which 
displays a genome-wide significant association of a 'large' effect 
size in anti-CCP positive RA with an OR above one (OR 1.45, 
95% CI 1.28 to 1.65, p=l. 33x10-^), representing a risk factor for 
this disease. However, no significant association is detected in 
anti-CCP negative RA (OR 0.91, 95% CI 0.73 to 1.13, p=0.37). 
The effect size ratio is highly significant (OR 1.60, 95% CI 1.27 
to2.02, p=7.36xl0-5). 

DISCUSSION 

With a number of anti-CCP-negative RA patients ranging from 
808 to 1918 per locus, our study represents the largest genetic 
study on anti-CCP negative RA to date. The 31 confirmed RA 
susceptibility loci identified so far have been primarily estab- 
lished in anti-CCP positive RA, and the meta-analysis by Stahl 
etaP-'^ included only seropositive RA. The effect size and strength 
of association of the loci identified in the meta-analysis by Stahl 
et al and in the anti-CCP positive RA subset of this study are 
very consistent. Slight differences might be related to the differ- 
ence of power for certain loci and to the use, in some instances, 
of several different proxies for one locus. It should be noted 
that we defined the minor allele frequency according to the 
frequency in the total population, including cases and controls. 
Since the minor allele frequency is close to 50% for SPREDZ in 
the meta-analysis by Stahl, where the minor allele frequency is 
based on controls only, the G allele is the minor allele, while it is 
the major allele in this study. This results in an inversion of the 
OR for association with anti-CCP positive RA. All loci tested 
here are associated with anti-CCP positive RA with p<0.05, and 
five loci reach genome-wide significance (<5E-08). 

Interestingly, we show a strong and highly significant asso- 
ciation of the SE with anti-CCP negative RA; the association 
remains present after stratification for rheumatoid factor (online 
supplementary table S2), which is consistent with the findings 
of a large study performed in patients from the same genetic 
background.'^ However, several smaller studies, performed in 
populations of different origins, did not detect significant associ- 
ation, likely due to a lack of power or ethnic differences. A large 
European meta-analysis, taking only 2-digit typing into account 
across several different European populations, concluded that 
there were weak, but no robust associations between any HLA- 
DRJBl alleles and anti-CCP negative RA.'^ The number of anti- 
CCP negative RA patients tested in our study is much higher 



Table 3 Schematic classification of RA susceptibility loci into three 
categories depending on their association pattern in anti-CCP positive 
and negative RA 



Category 

1 

2 
3 

Not classifiable 



Associations 

Both CCP positive and 
negative RA, stronger in CCP 
positive RA 

Both CCP positive and negative 
RA, equally strong in both 
CCP positive RA only, 
significant difference between 
CCP positive and negative RA 

CCP positive RA only, but no 
significant difference between 
CCP positive and negative RA 



Locus Name 

HLA-DRBl SE, PTPN22, 
TNFAIP3 locus 1 



AI\IKRD55, BLK, C5orf30, STAT4 

AFF3 locus 1, CCR6, CCL21 locus 
2, IL2RA locus 2, CD28, CD40, 
PXK, REL, RBPJ, TNFRSF14, 
TNFAIP3 locus 3 
All others 



than that in the meta-analysis or any other study, and confirms 
that the SE is associated with anti-CCP negative RA in UK 
Caucasians. 

Compared with the results in anti-CCP positive RA, where 
all loci show evidence for association, only seven loci achieve 
a p value below 0.05 in the anti-CCP negative subgroup. 
However, this number is higher than reported to date, and 
represents a statistically significant overlap (p=0.0018). Among 
the loci tested here, 32 non-HLA loci have been convincingly 
shown to be independently associated in anti-CCP positive RA. 
Therefore, without assuming any prior probability for associa- 
tion, the Bonferroni-corrected significance threshold would be 
1.6E-03. Three non-HLA loci would remain significantly associ- 
ated with anti-CCP negative RA after the Bonferroni correction: 
ANKRD55/IL6ST, BLK and PTPNZ2. 

Interestingly, both BLK and STAT4 are confirmed systemic 
lupus erythematosus (SLE) susceptibility loci and it would be 
interesting to investigate the overlap of other SLE loci with anti- 
CCP negative disease. Our findings confirm those of the previ- 
ous reports showing association in both anti-CCP positive and 
negative RA subgroups at the STAT4 locus. IRF$, which has been 
found to be associated with anti-CCP negative RA by previous 
investigators, is not associated with anti-CCP negative RA in 
this study. However, it should be noted that the current study 
only had 19.6% power to detect an association of an effect size 
of 1.12 at the 0.05 significance level. In addition, a differential 
association of this marker with RA and its different serotypes 
in different ethnic groups, or gene-environment interactions, 
might further account for the lack of association seen here in 
anti-CCP negative patients. 

The effect size of BLK and ANKRD55/IL6ST is larger in anti- 
CCP negative RA than in anti-CCP positive RA, although this 
difference does not reach statistical significance. All other loci 
show a higher effect size in anti-CCP positive RA. This can be 
explained by the fact that loci tested here have been primarily 
identified in cohorts containing mainly or exclusively anti-CCP 
positive patients. 



\ 1 



^1 



..I 



TNFAIP3 locus 3 



CCP, cyclic citrullinated peptide; RA, rheumatoid arthritis. 



Figure 1 OR and 95% CI for different single nucleotide polymorphism 
association patterns. PTPN22: significant association in both serotypes, 
significantly different. C5orf30: significant association in both 
serotypes with the same effect size. CCR6: only associated in anticyclic 
citrullinated peptide (anti-CCP) positive rheumatoid arthritis (RA). BLK: 
associated in both serotypes, effect size slightly, but not significantly, 
larger in anti-CCP negative RA. TNFAIP3 locus 3: only associated in 
anti-CCP positive RA, highly significant difference. PTPN22, C5orf30, 
CCRBare prototypic examples illustrating the three categories presented 
in table 3. 
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The genetic difference between the two RA serotypes can be 
particularly illustrated by the association pattern of the intronic 
locus at TNFAIP3: it is strongly associated with anti-CCP posi- 
tive RA, but not with anti-CCP negative RA, despite a detection 
power of over 80%; the difference in OR is statistically signifi- 
cant. The intronic locus at TNFAIP3 therefore represents a risk 
factor for anti-CCP positive RA, but not for anti-CCP negative 
RA in this analysis. 

The present study shows a genetic contrast between anti- 
CCP positive and anti-CCP negative RA and allows the classifi- 
cation of known RA susceptibility SNPs in different categories. 
The first category comprises markers associated with both sub- 
sets, but their effect size is significantly larger in anti-CCP posi- 
tive RA. Although anti-CCP negative RA could not be clearly 
divided into distinct clinical subphenotypes in a recent study in 
the Netherlands,^^ it might still comprise several genetically and 
serologically different subsets, based for example on the pres- 
ence of ACPA, other than anti-CCP antibodies.^' The second 
category contains SNPs, similarly associated in both anti-CCP 
positive and negative RA, with the same effect size. The third 
category comprises SNPs associated with anti-CCP positive 
RA, but not with anti-CCP negative RA and the effect size ratio 
is statistically significant. A lack of association of some mark- 
ers, while others are associated with disease irrespective of the 
serological status, suggests that RA susceptibility markers might 
cluster to different molecular pathways, some associated with 
autoantibody production, others not. 

In summary, among 33 independent genetic loci tested in 
this study, 18 could be classified into three different catego- 
ries, according to their association pattern, while 15 could not, 
mainly due to lack of power. Seven markers show an associa- 
tion with anti-CCP negative RA, while 1 1 others are unlikely 
to be associated. The use of a multinomial logistic regression 
analysis leading to three p values, as described here (anti-CCP 
positive RA, anti-CCP negative RA, effect size ratio), repre- 
sents a straightforward method to classify markers into three 
meaningful categories or to exclude them from classification, if 
only one p value out of three is significant. This latter situation 
occurs mainly when the power is low. The accuracy of classifi- 
cation depends on the definition of the significance threshold. If 
a Bonferroni corrected p value is used for classification, markers 
with p values between 0.05 and 1.6E-03 would change category. 
Markers with highly significant associations and effect size ratio 
like PTPN22 can be considered as accurately categorised. Future 
studies might identify more categories; for example, for makers 
associated exclusively with anti-CCP negative RA, or displaying 
a larger effect size in anti-CCP negative RA than in anti-CCP 
positive RA. 

Despite this being the largest sample of anti-CCP negative 
cases studied to date, the main limitation remains lack of power 
for many markers. This is particularly pertinent to the loci, 
which could not be classified into one of the three categories 
presented here. Larger sample sizes will be required to explore 
these loci more fully. Six thousand five hundred anti-CCP nega- 
tive patients and 1 1 000 controls would be required to detect an 
effect size of 1.10 with a power of 80% for a marker present at a 
MAP of 28% in controls (average MAP of RA susceptibility loci 
in controls reported in the study by Stahl et al^^). 

The low number of anti-CCP positive RA susceptibility loci 
associated with anti-CCP negative RA highlights the need 
for a well-powered GWAS for the discovery of yet unknown 
anti-CCP negative RA specific loci. The current study presents 
genetic differences and similarities between anti-CCP positive 
and anti-CCP negative RA. Although the two disease serotypes 



show significant differences in disease course and severity, anti- 
CCP antibodies are currently not used to guide treatment deci- 
sions in clinical practice. However, the results presented here 
highlight the need for genetic analyses of susceptibility, sever- 
ity and treatment response to consider the two serotypes both 
separately and together for future investigations. 
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